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Abstract: We analyze the spectrum of massive scalar bound states in the background of 
extremal Kerr black holes, focusing on modes in the superradiant regime, which grow ex- 
ponentially in time and quickly deplete the black hole's mass and spin. Previous analytical 
estimates for the growth rate of this instability were limited to the jiM <C 1 and fiM S> 1 
regimes, where /i and M denote the scalar field and black hole masses, respectively. In this 
work, we discuss an analytical method to compute the superradiant spectrum for generic 
values of these parameters, namely in the phenomenologically interesting regime ~ 1. 
To do this, we solve the radial mode equation in two overlapping regions and match the 
solutions in their common domain of validity. We show that matching the functional forms 
of these functions involves approximations that are not valid for the whole range of scalar 
masses, exhibiting unphysical poles that produce a large enhancement of the growth rate. 
Alternatively, we match the functions at a single point and show that, despite the un- 
certainty in the choice of the match point, this method eliminates the spurious poles and 
agrees with previous numerical computations of the spectrum using a continued-fraction 
method. 
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1. Introduction 

The issue of black hole stability plays a central role in General Relativity and is crucial 
for our understanding of the rich high-energy phenomenology associated with compact 
astrophysical systems. The issue of stability of the Schwarzschild geometry was first ad- 
dressed in the seminal paper by Regge and Wheeler [1] and later by Vishveshwara [2] and 
Wald [3] , showing that small perturbations of this spherically symmetric space were expo- 
nentially damped in time. It has been known for a few decades that this is not, however, 
the case for rotating Kerr black holes [4], which suffer from the well-known superradiant 
instability [5-13]. 

This instability is a direct consequence of the presence of an ergoregion surrounding the 
outer horizon of the black hole, where ordinary causal matter experiences inertial frame- 
dragging and cannot remain at rest with respect to an asymptotic observer. Within this 
region, the Killing vector associated with time translations becomes space-like, so that 
the associated conserved energy may become negative for some observers. Classically, this 
gives rise to the so-called Penrose process [14,15], by which one may extract energy and 
angular momentum from the spinning black hole. One can consider a projectile directed 
towards the black hole along a time- like trajectory but aimed in such a way that it misses 
the outer horizon. One could then conceive a timer mechanism that sets off a fragmentation 
process that breaks the projectile into two parts while it is inside the ergoregion. One of 
the fragments may then carry a negative energy and fall into the black hole, while the 
remaining fragment escapes to infinity carrying an energy larger than that of the original 
fragment. One can also show that the same reasoning holds for the angular momentum of 
both fragments. 

Superradiance is a well-known phenomenon in several quantum [16, 17] and also clas- 
sical systems [18] and was first described for the case of spinning black holes by Zel- 
dovich [5,6]. If one considers the scattering of wave modes of the form e~*'^*+*'"''^ , where (/> 
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denotes the Boyer-Lindquist azimuthal angle in the Kerr geometry [20] , one can show that 
for co-rotating (m > 0) incident waves satisfying the condition oj < m^l, where = •2°, 2 
is the angular velocity of the black hole at the outer horizon r+, the scattered wave is 
amplified. The total energy and angular momentum radiated to infinity are hence larger 
than that of the incident wave, in a similar fashion to what happens in the Penrose process 
described above. Press and Teukolsky proposed that, by surrounding the black hole with 
a mirror, one could use multiple wave-scattering to extract an enormous amount of energy 
from the black hole, a mechanism they suggestively named the black hole bomb [21] and 
that was further analyzed in [19]. 

Practical implementations of such a mechanism seem difficult to conceive, although 
it has been proposed in [22, 23] that the inner boundary of the accretion disks which are 
expected to surround astrophysical black holes may effectively act as a mirror and pro- 
duce enough energy via superradiant scattering of magnetosonic plasma waves to feed the 
mysterious and highly-energetic gamma-ray bursts. A natural mirror is, however, provided 
for fields with a non-zero mass [7, 8] . It is well-known that stable orbits around a black 
hole are an exclusive feature of particles with a non- vanishing mass, so that massive fields 
exhibit a set of bound states in the Kerr background corresponding to wave-packets mov- 
ing along such orbits. These wave-packets will always have a non- vanishing tail inside the 
black hole's ergoregion, leading to a continuous amplification of bound states satisfying the 
superradiant condition. Such states have a complex frequency with a small positive imag- 
inary part, so that their occupation number grows exponentially in time. Pauli blocking 
prevents this amplification for fermionic states [24], so that superradiant bound states are 
an exclusive feature of bosonic fields. 

Several attempts have been made in the literature to analytically determine the growth 
rate of superradiant bound states for massive scalar fields. The difficulty in solving the 

massive Klein-Gordon equation in the Kerr background has limited the analysis to the 
regimes Af/x ^ 1 [8] and Mfi <^ 1 [9, 10], where M and fi denote the masses of the black 
hole and the field, respectively, in geonietrized units such that G = c = h = 1. In both 
cases the obtained growth rates were shown to be quite small, with the imaginary part 
of the frequency cojM scaling as e"^-^^'*^ and {fiM)^ in the large and small mass limits, 
respectively. These studies suggest as well that the growth rate is maximal for extremal 
holes with a = M and for P- waves with I = m = 1. The spectrum of massive scalar bound 
states has also been studied numerically in [12] and later in [13], using a continued- fraction 
method which confirmed these expectations. The latter analysis revealed a maximum 
growth rate in the I = m = 1 case given by ujM ~ 1.5 x 10^'', for /xM ~ 0.42 and 
a = 0.999M, in agreement with the former, where a similar value was observed for a = 0.994 
and nM = 0.45. 

Hence, numerical results suggest that the superradiant instability is more pronounced 
in the regime /xM ~ 1, which has proven hard to analyze with analytical methods, mo- 
tivating a recent analytical study of this phenomenon [25]. By separating the exterior of 
the black hole into two overlapping regions where the field equations are exactly solvable 
and matching the obtained solutions within their common domain of validity, the authors 
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obtained a defining equation for the superradiant bound state spectrum valid for ~ 1. 
The novelty of their method resides in the use of the black hole's angular momentum rather 
than the particle's mass to define the two regions, thus eliminating the need for large and 
small mass approximations in order to solve the relevant equations. The results in [25] 
exhibit, however, a huge discrepancy with respect to the numerical analysis in [12,13], 
with a maximum growth rate of ujM 2± 1.7 x 10~^, four orders of magnitude larger than 
the one mentioned above. 

Understanding this discrepancy is crucial for the study of the phenomenological aspects 
of scalar superradiant bound states in the Kerr background, in particular as this process 
may lead to the significant spin down of rapidly-rotating black holes. Astrophysical black 
hole candidates have been primarily found in X-ray binaries and Active Galactic Nuclei 
(AGN), with masses in the range 3Mq — 30Mq for stellar mass black holes and 10^ Mq — 
10^ Mq for supermassive black holes, where Mq denotes the solar mass [26]. The lightest 
known scalar particle is the neutral pion, with n = 134.98 MeV and a mean lifetime 
c± 8.4 X 10"-'^^ sec [27], so that it may only significantly affect small black holes with 
a mass of about 10^^ kg, given the requirement fiM ~ 1. Hence, pionic bound states are 
not expected to produce observable effects on astrophysical systems, although they may 
influence the dynamics of small primordial black holes postulated to be produced during 
inflation [8]. This statement depends crucially on whether the instability can develop 
significantly before the pions decay. While a large amplification of pionic bound states can 
be attained using the analytical results of [25], where the time scale for the instability is 
T ~ 1.3 X 10~^^ sec, the numerical results of [12, 13] yield r ~ 1.5 x 10~^^ sec, in which 
case superradiant pion emission would be very ineffective. 

Recently, the authors of [28] considered the many pseudoscalar moduli fields with 
axion-like properties which generically appear in string theory compactifications. These 
moduli arise as Kaluza-Klein zero-modes of the Neveu-Schwarz and Ramond-Ramond 
closed string antisymmetric form fields and, if massless at tree level, may acquire masses 
from non-perturbative effects such as string instantons. The authors then argued that, if 
realistic compactifications are to include the QCD axion as a solution to the strong CP 
problem, string instanton contributions to its potential must be suppressed with respect 
to non-perturbative QCD contributions by at least a factor of 10^'^. This suggests that 
many of the axion-like moduli will be lighter than the QCD axion, generating a plethora of 
ultra-light scalar states they denoted as the "string axiverse" , with masses possibly down 
to the Hubble scale, Hq ~ 10~^^ eV. Many of these axions may thus lie in the correct mass 
range to significantly spin down astrophysical black holes, possibly leading to several inter- 
esting observational effects such as axion-photon conversion, through the coupling to the 
surrounding magnetic fields, and emission of gravitational waves, via quantum mechanical 
transitions between superradiant and non-superradiant states. 

Spinning astrophysical black holes may thus provide a simple probe of string theory 
compactifications. In particular, black holes with mass M ^ for a given axion should 
exhibit a very low angular momentum, so that looking for gaps in the mass-spin Regge 
plot for astrophysical black holes may provide an effective method to search for ultra-light 
states [28]. However, the shape of these gaps will be extremely sensitive to the variation 
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of the superradiant growth rate with the scalar mass. Moreover, black holes may also gain 
angular momentum from the surrounding accretion disks, a process that may compete with 
the superradiant spin down if the instability does not develop sufficiently fast. 

The rich phenomenology associated with the superradiant instability thus requires a 
more detailed study of its properties and an analysis of the above mentioned discrepancy. 
With this goal in mind, in this paper we compute the spectrum of superradiant massive 
scalar bound states in an extremal Kerr background, for which the growth of the instability 
is expected to be maximal. We consider a functional matching procedure analogous to the 
one used in [25] and show that, unfortunately, it is not applicable for the whole scalar 
mass range, mainly due to the approximations involved. Alternatively, we consider a point 
matching procedure and show that, within the uncertainty associated with the choice of 
the match point, this method agrees with the numerical results of [13] for a large range of 
masses, in particular in the vicinity of the above mentioned discrepancy. 

This work is organized as follows. In the next section, we analyze the main properties 
of massive Klein-Gordon fields in a maximally-rotating black hole. Sections III and IV 
are devoted, respectively, to the analyses of the functional and point matching procedures. 
We conclude, in Section V, with a discussion of the main advantages and disadvantages of 
both techniques and some prospects for future work. 

2. The Klein-Gordon equation in an extremal Kerr background 

We start by considering the Klein-Gordon equation in the Kerr background, which for a 
scalar field <I> of mass /i can be written as: 

(V^V^ - = . (2.1) 

In terms of the Boyer-Lindquist coordinates {t,r,6,(f)) [20], the scalar field admits a mode 
expansion of the form: 

•^Ut, r, e, <P) = e-"^'RUr)SUe)e'""^ , (2.2) 

where u = (jJr + iujj is the complex frequency of the mode characterized by the indices 
/ = 0, 1, . . . and —l<m<l. The latter correspond to the usual angular momentum 
and azimuthal projection indices for spherical harmonics in the Schwarzschild limit a — )• 
0. In the general case, they label the solutions of the spheroidal harmonic equation [29] 
corresponding to the angular part of Eq. (|2.1| ): 
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where q = -yZ/x^ — and Aim is the angular eigenvalue, corresponding to the separation 
constant relating the radial and angular parts of the Klein-Gordon equation. The regularity 
of the solution of Eq. ( ^ ) at the poles 9 = 0,tt determines a discrete set of angular 
eigenvalues, which for a^q^ <C 1 are given by an expansion of the form: 

oo 

Aim = l{l + l)+Y, Cklma^'^q^'' . (2.4) 
k=l 
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We are interested in studying bound states for which uj is only shghtly smaller than ^u, i.e. 
q 1, so that it will be sufficient to consider the first few terms in this expansion, which 
can be found in [29] . The radial Teukolsky equation is given by [30] : 



dr \ dr 



2 2 2 2/1 

— a (jj + zmacu — ji r — Aim 



Rim = 0, (2.5) 



where A = r^ + a^ — 2Mr and K = {r'^ + a'^)uj — am. The zeros of A at r± = M^s/M"^ — 
correspond to the outer and inner horizons of the Kerr black hole. 

For extremal Kerr black holes, with a = r+ = r_ = M, it is useful to introduce the 
variable: 

and the re-scaled function ^'/^(x) = xi?;^(x), which satisfies the Schrodinger-like equation: 

^^+[u;^-V{u;,x)]^lm = 0, (2.7) 



where 



V{u:,x) = -—^ ^ + ^2 (2-^ 



and we have defined: 



2^2 _ 2 

u = , w = uj — mQ , (2.9) 



P{/3-l) = A- 7uj^ + , (2.10) 

where all quantities are measured in units of M, Q = 1/2 in the extremal case and we have 
omitted the indices l,m for simplicity. Without loss of generality, we use the following 
solution for the coefficient (3: 

/3 = ^ (^1 + yTT4J^^2&j2^V j ~ i + 1 , (2.11) 

where we have taken the limit w ~ ^ 1 in the last expression, explicitly showing that 
f3 labels the angular momentum of the solution. It is a particular feature of the extremal 
case that one can readily write the radial equation in this form without using the tortoise 
coordinate, defined by [8]: 

dy r'^ + ,2 , , 

1 = ^ ^ y = x + 21og.--. (2.12) 

Despite its simple form, the potential determining the dynamics of the radial function 
depends on the mode frequency, making this problem hard to solve exactly. To better 
understand the form of this potential, we use the solution for the real part of the mode 
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Figure 1: Potential for I = m = 1 and fiM = 0.4 (solid line). Region I corresponds to a; < x^ax — 
(2g)-i, while region II is defined by x > x„iin — — 4n7. The point a;, = ^ —Iw jq corresponds to the 
trial match point discussed in Section IV. The dashed and dotted curves give approximate forms 
of the potential valid in regions I and II, respectively. 



frequency that we will obtain later on in Eq. (|3.5| ). This is illustrated in Figure 1 for 
/ = m = 1 and [iM = 0.4. 

The potential can be divided into (i) a sharp potential close to the horizon, (ii) a finite 
centrifugal barrier and (iii) a potential well. In terms of the tortoise coordinate, the sharp 
part of the potential tends to a constant value 1m£lL>j — rr?^ close to the horizon, which 
is negative for superradiant modes [28]. The definitions of regions I and II, as well as the 
particular values of x shown in Figure 1, will become clear in the next sections, but they 
roughly separate the interior and the exterior of the ergoregion. 

The physical picture behind the superradiant amplification of the modes has a simple 
interpretation in terms of this potential. The Hydrogen-like bound states of the potential 
well may "anti-tunnel" [28] into the ergoregion and, while a fraction may penetrate the black 
hole's horizon, the remaining part will be reflected with a larger amplitude for a; < m£l. 
This energy remains localized in the potential well rather than being radiated to infinity, 
as the wave packets are bound to the black hole, effectively inducing multiple scatterings as 
if a mirror was placed around the horizon. This results in a small positive imaginary part 
of the frequency, making the occupation number of the bound states grow exponentially 
according to Eq. (|2.2| ). This effect should be smaller for higher multipoles, for which the 
barrier height increases and the "tunneling amplitude" is suppressed. 

In order to determine the spectrum of bound states, one needs to impose boundary 
conditions. These correspond to an "ingoing wave" at the horizon and an exponentially 
decaying mode at infinity [8, 10], which in terms of the tortoise coordinate defined in Eq. 



(2.12) can be written as: 



^Uy)^\~'"~^^'' ^ . (2.13) 
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Recall that y ~ — § and y ~ x close to the horizon and at infinity, respectively. One 
can rewrite the radial equation in a more convenient form by considering the function 
Fim{x) = x~^^im{x), which satisfies 



n cP ^ d 4tu^ Scow 9 9 

x^—r + 2/3x— + + + 2gzyx - q^x^ 

ax^ ax x'^ x 



. (2.14) 



This cannot be solved exactly but admits approximate solutions in two different regions. 
These may be matched in a common domain of validity to yield the spectrum of bound 
states with the boundary conditions in Eq. ( |2.13| ) , which will be the topic of the remainder 
of this paper. 

3. Functional matching 



As mentioned in the previous section, Eq. ( 2.14 ) can be simplified in two distinct regions 



by noting that the first two terms in the potential may be neglected for — 4tu/x <^ 1, where 
tz7 < for superradiant modes, while the last two terms may be neglected for 2qx <^ 1. 
This defines the regions I and II depicted in Figure 1 such that x < Xmax = (2(7)"^ and 
X > Xrnin = — 4tu, respectively. In Figure 1 we have also plotted the approximate forms of 
the potential obtained by neglecting the relevant terms in each region, denoted by Vi{x) 
and Vii{x). 

If one is able to determine solutions of the radial equation in these regions, one may 
determine the spectrum by matching the obtained functions in the overlap region, x^in ^ 
X <^ Xmax- The existence of such a region requires —8wq < 1, a condition that we cannot a 
priori ensure to be satisfied without computing the spectrum. We will nevertheless assume 
that this is the case and use it as a consistency check of the matching methods. 

The equations in both regions can be reduced to the confiuent hypergeometric equa- 
tion, although in region I this requires a change of variable to y = —2/x. Note that this 
variable corresponds to the near-horizon form of the tortoise coordinate, so that we use 
the same notation for both to emphasize this fact. Also, as y ~ x for large x, we can 
say that in both regions the radial equation reduces to the confiuent hypergeometric equa- 
tion in the corresponding limit of the tortoise coordinate. This reduction requires writing 
^(y) = e"^G(x) with an appropriate choice of the complex parameter a in each case. After 



imposing the boundary conditions Eq. ( |2.15 ), we obtain the following solutions: 

4:iTu' 



Fi{x) = Aie—U\^l-l3-2iuj,2-2l3, 
Fn{x) = Ane-'^'^Ui^ - u, 2/3, 2qx) , (3.1) 

where Ajji are normalization constants and U{a, b, z) is the confiuent hypergeometric 
function which is regular as \z\ — )• oo. In the overlap region, we may take the \z\ — )• limit 
of the latter, given by [29]: 



U{a,b,z] 



sin(7r6) 



1 z^- 



r[l + a-6]r[6] r[a]r[2-6] 



(3.2) 
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This yields for the corresponding ^ functions: 



^^lix) = -Ai ^ 



sin(27r/3) 



■^n{x) = An 



TT 



sin(27r/3) 



r[/3 - 2ia;]r[2 - 2/3] r[l - /3 - 2ia;]r 



r[i - /3 - z.]r[2/3] r[/3-z.]r[2-2/3] 



(3.3) 



Hence, both functions have a common functional form in the overlap region, and one may 
match the coefficients of the x^ and x^~^ terms to obtain the following condition: 

r[/3-z.] ^'"^"^^ r[l-/3-2ia;]V r[2/3] ' ^"^'^^ 

For the matching to be possible, —Swq < 1 and, as in the small mass limit 2/3 — 1 ~ 2/ + 1 > 
0, we conclude that the RHS of this condition is suppressed, so that the LHS should vanish 
to leading order, which implies: 

fi-.^-n^ ...,(^,-^jEMl_y ,3.5) 

where n is a non- negative integer. This reproduces the Hydrogen-like spectrum obtained 
in previous works for fxM <C 1 [10, 13], although significant deviations are observed for 
fj,M ~ 1, where corrections to the real part of the frequency arise both from the polynomial 



equation (3 — u = —n and from the finiteness of the RHS of Eq. (3^). This expression 
confirms nevertheless our physical intuition about the Hydrogen-like nature of the massive 
bound states in the black hole's potential well. 

We need to go beyond this leading approximation to obtain the imaginary part of the 
spectrum and with this purpose we have used Mathematica to numerically obtain the roots 
of Eq. ( |3.4D . The FindRoot command can be used for this purpose, requiring an initial 
value for which we considered the relevant root oi (3 — v = —n in each case. In Figure 2, 
we plot the results for I = m = 1 and n = 0, for which we expect the maximum growth 
rate. 

In Figure 2(a), one observes a smooth increase of the growth rate with ^u, with a sharp 
decrease close to the endpoint of the spectrum, where the superradiant condition is no 
longer satisfied and which for <^ 1 lies close to m/2. Although this is in agreement 
with the numerical results of [12,13], one also observes two peak-like structures separated 
by a sharp decrease in uJiM. One of these "peaks" was also observed in [25] and is behind 
the four orders of magnitude discrepancy between the two types of analysis. Figure 2(b) 
depicts a detail of this feature and should be compared with Figure 1 in [25]. One should 
note, however, that the matching condition Eq. ( |3.4| ) involves several gamma functions, 
some of which may develop poles for some values of /i. In particular, r[2 — 2/3] has poles 
for /3 = 1 + p/2, where p is a non- negative integer. Using the lowest order result in Eq. 
( |3.5| ) for w/j, one obtains poles in this function for: 

. (^'^^);-(^^')- , (3.0) 
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Figure 2: Results obtained for ttie imaginary part of the frequency obtained using the functional 
matching approach for I — m — I and 71 = 0, illustrating (a) the smooth part and the two peaks 
and (b) a detail of the second peak. 



which gives /u ~ 0.46 for I = p = 1, in the region where the peak- like features in the 
superradiant spectrum are observed. The shape of this pole is illustrated in Figure 3, 
showing that it occurs exactly at the value of /uM where the sharp decrease in the growth 
rate is observed. This suggests that the peak-like structures are most likely unphysical and 



simply a consequence of the approximations involved in deriving Eq. (3.4). 



In fact, the limit taken in Eq. (3.2) is not valid close to these poles. To see this, we 
may rewrite the solution in region I as: 



Fj{x) = Aje' 



2/3-1 



U[ /3 -2iuj,2(3,- 



where we used that [29]: 



[/(a, 6, z) = z^-^U{l + a-b,2-b,z) 



(3.7) 



(3.8) 



In this form, we have b = p + 2 for both solutions, p G Zq, so that the correct jz] — )■ 
limit of the confluent hyper geometric function is given by [29]: 



, , (-l)''[logz + ^(a)+7-^(&)] , r[b-l] 
U [a, b, z) ~ ^r.T^r T—r-, H ^rr-; — z 



r[6]r[a-6 + l] 



T[a] 



(3.9) 
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Figure 3: Real part of the function (r[2 — 2/3])^ appearing in the matching condition Eq. (3.4), 
illustrating the pole at the value of fxM where the sharp decrease in the growth rate of the instability 
is obtained using the functional matching method. 



where 'il^{x) denotes the Digamma function and 7 the Euler-Mascheroni constant. Using 
this expression, the solutions in the two regions become: 



ipn{x) 



Ai 
An 



{iiw)^l^-^[\og{-iiw/x) + ij{/3 - 2iu) + 7 - ^(2/3)] 

T[2/3]T[1 - /3 - 2iuj] ^ ' T[/3 - 2iuj] 



-1)^/^ [ log(2gx) + + 11^(2(3)] 

r[2/3]r[i -f3-u] 



r[/3 - z.] 



(3.10) 



Hence, the functions cannot be matched in this the log X corrections are associated 

with distinct powers of x in each solution. Such logarithmic terms are generically significant 
in the vicinity of points for which /? = 1 +p/2, so that the functional matching procedure 
outlined in this section should only be used away from the above mentioned poles. 

In particular, one cannot a priori neglect the effects of the logarithmic corrections near 
the peaks illustrated in Figure 2, even though all gamma functions are finite in this case, 
due to their obvious proximity to the pole. Furthermore, the discrepancy with respect to 
the numerical results of [12, 13] suggests that these corrections are indeed significant. The 
number of such poles increases as we consider higher multipoles, making it hard to trust 
this technique in general. This is a disappointing result, as it precludes the use of the 
functional matching method to cover the whole range of scalar masses. 



4. Point matching 

As we have seen in the previous section, the log-corrections preclude the matching of 
the functions in the overlap between regions I and II. This is not, however, the only 
possible technique one can use to obtain the spectrum, and a simple alternative consists in 
matching the functions and their first derivatives at a point. This technique is widely used 



-10- 



in problems in which the potential has two clearly separated regions where the solutions 
take distinct forms. In this case, however, we have two overlapping regions where one can 
find approximate solutions to the radial equation. Nevertheless, choosing a match point 
in the overlap of regions I and II, we may ensure that both functions are sufficiently close 
to the exact solution and thus obtain a good approximation to the spectrum. One cannot, 
however, hope to obtain an accurate prediction for the imaginary part of the frequency, as 
this would require an exact prescription for the choice of the match point, which we do not 
have a priori. 

Matching the functions in Eq. ( |3.lD at a generic point x, one obtains the following 
condition: 



2 

qx 
2iw 



l + 2(/3 



C/(l - /3 - 2iuj, 2 - 2/3, -4iw/x) 
, [7(^-1/ + 1,2/3 + l,2gx) 



C/(/3 - u, 2f3,2qx) 



(4.1) 



where we have used that U'{a, b, z) = —aU{a + 1, 6 + 1, z) [29]. A good trial match point 
is X* = \J —2w jq^ as this corresponds to the value of x for which the absolute values 
of the arguments of the confluent hyper geometric functions in both regions are the same, 
simultaneously giving the geometric mean of the values of Xmin and Xmax defined in Section 
II, which give the boundaries of the two regions. This ensures that the matching condition 
Xmin < X < Xmax is Satisfied if such an overlap region can be found. We then obtain for 
this particular choice: 

l + 2(l-/i-2i.)'^i^-^^'"^^-^^■'-'•!=.fl + 2(g-.)^'^-^' + '^^^+^-'•' 



t/(l -^-2ia;, 2-2/3, iz,) V C/(/3 - i/, 2/3, z,) 

(4.2) 

with = ^ySwq < 1 for an overlap region to exist. We have determined the roots of 
this condition numerically for I = m = 1, 2 and 3, which give the largest growth rates, 
via the same procedure used for the functional matching method in the previous section. 
These results are plotted in Figure 4, along with the corresponding solutions obtained 
using functional matching and the numerical continued-fraction method in [13] for the case 
a = 0.999M, which corresponds to the largest black hole spin analyzed in the latter work. 
As mentioned above, there is no absolute prescription for the choice of the match point, 
so that we also plot in this figure the results obtained by choosing different points within 
the overlap region, Xmin < x < Xmax, giving an effective uncertainty for the match point 
method. 

As one may observe in the three plots in Figure 4, the match point procedure gives a 
smooth evolution of ujM as a function of fiM and does not exhibit any peak-like features as 
those obtained using the functional matching. For the first three multipoles, these results 
are in good agreement with the numerically obtained values for a = 0.999M, although they 
differ significantly close to the endpoint of the superradiant spectrum. In fact, the endpoint 
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Figure 4: Results obtained for ttie imaginary part of the mode frequency using the point matching 
method for Z = r7i = 1, 2 and 3, with n = 0, along with the uncertainty associated with the choice 
of match point. Also shown are the results obtained using the functional matching method and the 
numerical solutions of [13] for a = 0.999Af. 



corresponds to = 0, which is extremely sensitive to the black hole's spin, extending the 
superradiant spectrum to larger /i as one approaches extremality. 
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For I = m = 1, both the numerical and functional matching results (away from the 
pole) lie within the shaded region and are in good agreement with the values obtained from 
point matching at x*, showing that in this case the latter is extremely close to the optimal 
match point. The upper (lower) limit of the shaded region corresponds to matching at 2^rnm 
{xmax), mainly due to the fact that the Hydrogen-like function 'tljjj{x) has a larger overlap 
with the ergoregion if the matching is implemented closer to the horizon. The maximum 
growth rate for matching at x* occurs for fj,M ~ 0.454, for which u!jM ~ 1.49 x 10"'^, in 
good agreement with the numerical results of [12,13], despite the difference in the black 
hole's angular momentum and the uncertainty in the choice of the match point. These 
results clearly show that the four orders of magnitude discrepancy between the results 
of [25] and [12, 13] are unphysical. 

For I = m = 2 and 3, Figure 4 shows only a small range of /iM, mainly because for large 
masses the number of poles in the functional matching condition becomes quite large, which 
precludes any decent comparison with the point matching technique. For the mass values 
shown in Figures 4(b) and 4(c), one observes a significantly larger discrepancy between the 
two analytical techniques away from the poles with respect to the I = m = 1 case, with 
the point matching results lying much closer to the corresponding numerical curves. This 
suggests that the approximations involved in deriving the functional matching condition 
fail for these multipoles even for /iM away from the expected poles, so that this method 
does not give reliable results in this case. 

To get a better comparison between the point matching and the numerical methods, 
we plot in Figure 5 the results for the first three multipoles with / = m and n = 0, including 
values for matching at x*, the associated uncertainty for the whole matching region and 
also the numerical curves for a = 0.99M and a = 0.999M obtained in [13]. 
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Figure 5: Results for the imaginary part of the mode frequency obtained using the point matching 
technique for the first three multipoles with I — m and n — 0, including the corresponding numerical 
solutions of [13] for a = 0.99M and a = 0.999Af . 

As expected, the maximum growth rate decreases for larger values of / = m, as in 
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these cases the overlap of the radial function with the ergoregion becomes smaller. The 
agreement between the point-matching and the numerical curves is extremely good, apart 
from the vicinity of the endpoint where one observes a large discrepancy. In fact, these 
deviations are larger for the higher multipoles, which is expected as the sensitivity of the 
condition w = to the value of a is much greater for larger m. This can be readily seen by 
comparing the two numerical curves for a = 0.99M and a = 0.999M for each multipole. 

The spread in the values of ujjM due to the size of the matching region is larger in 
the neighborhood of the maximum growth rate for the first two multipoles. As discussed 
earlier, this can be seen as an effective uncertainty associated with the point matching 
technique, although one should take into account that both Vj{x) and Vji{x) are not good 
approximations to the radial potential close to the limits of the matching region. Also, a 
larger spread in the values of ujM also indicates a wider overlap region, ensuring that both 
ipi{x:t) and ^pn{x:^) are good approximations to the exact radial function at x*, which is 
confirmed by the good agreement with the results for a = 0.999M up to the value of 
giving the maximum growth rate. 

On the other hand, the results for I = m = 3 exhibit a more significant deviation from 
the numerical ones for small fiM. This is somewhat puzzling, as the uncertainty due to the 
choice of the match point is negligible in this region. To better investigate these results, 
we have used the spectrum computed at the match point x* to determine the boundaries 
of regions I and II for this multipole. The separation of the two regions for this particular 
case is illustrated in Figure 6. 
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Figure 6: Matching region for I = m = 3 and n = as a function of the scalar mass. The dotted 
and dashed curves correspond to Xmin — —4-07 and Xmax ~ (2'?)^^, the lower and upper limits of 
regions II and I, respectively. The solid line corresponds to the trial match point a;, ~ ^ —2v3 j q. 

As one can easily conclude from this figure, there is no overlap region for 0.8 < [iM < 
1.4, which explains the more significant discrepancy for I = m = 3 with respect to the 
numerical results in this part of parameter space. For such values of nM, Xmin is only 
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slightly larger than Xmaxi so that the spread in the corresponding values of ojjM is quite 
small. One has to take into account, however, that some part of this deviation is also 
due to the non-extremality of the numerical results, a fact that precludes a more accurate 
determination of the effects of the absence of a matching region. Nevertheless, we may safely 
conclude that, despite this difference, the point matching at gives a good approximation 
to the numerically obtained values of cojM. 

5. Conclusions and outlook 

In this work, we have analyzed the spectrum of superradiant bound states of a massive 
scalar field in the background of an extremal Kerr black hole. We have obtained a defining 
equation for this spectrum using two different techniques, both based on the assumption 
that one can separate the exterior of the black hole into two overlapping regions where 
the radial equation is exactly solvable. We have then computed the numerical roots of 
these conditions in order to determine the growth rate of the superradiant instability as a 
function of mass of the scalar field. 

The first technique, analogous to the one used in earlier works such as [10] and [25], 
is based on matching the functional forms of the two solutions in their common domain 
of validity. We have shown that, although this technique is valid for a large range of 
scalar masses, it exhibits anomalous points which correspond to poles in the matching 
condition and in the vicinity of which the approximations involved do not hold. For I = 
m = 1, one finds two peak-like structures in the vicinity of one of such points, which had 
been interpreted in [25] as a physical enhancement of the growth rate of the instability, 
disagreeing with earlier numerical computations [12, 13] by four orders of magnitude. 

As an alternative to the functional matching method, we propose matching the two 
functions and their derivatives at a single point within the overlap region. Choosing the 
match point at the geometric mean of the limits of this region, = ^/—2wjq, we obtain 
a good agreement with the numerical results of [13]. The non-extremality of the latter 
precludes a more accurate comparison of the results obtained using both techniques close 
to the endpoint of the superradiant spectrum at a; = m^l, but the agreement for small 
liM suggests that the results obtained for can be trusted for the whole range of scalar 
masses. In fact, one could use the numerical results for small masses to better calibrate the 
choice of the match point, although the results obtained for x* are sufficient to illustrate 
the main properties of the scalar superradiant spectrum. 

This method also reveals no anomalous points but only a smooth increase of the growth 
rate with the scalar mass, with the expected sharp decrease close to the endpoint. Despite 
the absence of a precise prescription for choosing the match point, this method has a 
significant advantage over the functional matching method, as it allows a semi-analytical 
study of the instability over a large range of masses. 

We have analyzed the first three multipoles with / = m, observing the expected de- 
crease in the maximum growth rate of the instability with increasing I. We have obtained 
a maximum value tJjM ~ 1.49 x 10~^ for /iM ~ 0.454 in the I = m = 1 case, in close 
agreement with the results of [12, 13]. 
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For pions, this corresponds to a time scale for the development of the instability of 
r = 1.43 X 10"^'' sec in a Kerr black hole with M = 8.65 x 10^^ kg. This is only about 
17% of the pion's mean lifetime, which suggests that, even for small primordial black holes, 
pions cannot extract sufficient energy and angular momentum from the black hole via 
the superradiant instability. The proposed (string) axion-black hole bomb [28] remains, 
however, a viable possibility for astrophysical black holes and it is worth investigating the 
rich phenomenology associated with it. 

Both the functional and point matching procedures depend crucially on the existence 
of an overlap between the near and far regions, which for this analysis corresponds to the 
condition —^wq <^ 1. This implies that the matching methods are valid either close to 
the endpoint of the superradiant spectrum, where w = uj — <^ 1, or for slightly bound 
states, for which q = \J p? — oS^ <^ 1. While for the lowest multipoles at least one of these 
conditions is satisfied for the whole range of scalar masses, in the case I = m = 3 there is 
a region of parameter space where no overlap between the two regions exists, effectively 
rendering any matching procedure inconsistent. This is mainly due to the fact that the 
lower limit of region II, Xmin = — 4ro, significantly increases for larger azimuthal number 
m. We cannot a priori quantify this increase, as it depends on the mode frequency, but 
one should expect this to be a generic feature of higher multipoles. This does not have, 
however, important phenomenological consequences, as the growth rate of the superradiant 
instability is significantly smaller for these multipoles. Also, despite the absence of an 
ovelap region, the difference between the point matching results and the numerical results 
of [12, 13] is small for a large range of scalar masses, so that the former may be sufficiently 
accurate for most phenomenological purposes. 

The existence of an overlap region for all the values of /xM considered for I = m = 
1 ensures the applicability of the point matching procedure in this case. On the other 
hand, the analysis of [25] is restricted to the regime u ~ mQ ~ /i, precluding the use of 
the matching procedure for small scalar masses. This difference arises from the distinct 
definitions of the regions I and II (or near and far regions, respectively) used in this work 
and in [25], with the extremality condition playing an important role in this case. 

It is worth mentioning that the results obtained in this work, as well as the numerical 
analysis in [12, 13], seem to be in confiict with those of [11], as pointed out in [31]. In this 
numerical analysis of the scalar superradiant instability, the authors estimated a growth 
rate ujM = 2 x 10^^ for fiM = 0.25, which is two orders of magnitude larger than the 
maximum value obtained in the present work. The numerical method used in [11] is, 
however, fundamentally different from the continued- fraction method employed in [12,13]. 
In particular, it computes the time evolution of a particular scalar mode in the Kerr 
background rather than determining the spectrum of superradiant bound states. Moreover, 
the selected mode has a frequency w = /i which, as discussed in the present and in earlier 
analyses, does not correspond to a scalar bound state. Although our results show that 
the frequency of bound states is very close to the scalar mass /i when the latter is small, 
as explicitly obtained in Eq. (|3.5| ), they also suggest that small deviations in the real 
part of the spectrum may produce large (unphysical) effects in the growth rate of the 
instability. Also, the above mentioned value is derived under the assumption that the 



-16- 



mode is growing exponentially, while the estimated e-folding time is in fact much larger 
than the numerically sampled interval. Thus, the claim in [31] for the existence of a regime 
where the superradiant growth rate is much larger than the results obtained in this work is, 
in our opinion, unjustified and involves a non-trivial extrapolation of the results obtained 
in [11]. 

In the comment to the analysis described in this work of [31], the authors of [25] 
explicitly show that all gamma functions involved in their (functional) matching condition 
are finite for values of fj,M in the region corresponding to one of the peaks obtained for 
I = m = 1. This is in agreement with the results obtained in this work, where r[2 — 2/3] 
exhibits a pole in between these peaks but is finite in the peak regions, as shown in Figure 
3. However, as discussed in Section 3, the existence of a pole in the vicinity of the peaks 
precludes the use of the form in Eq. ( |3.2D for the \z\ — )• limit of the hypergeometric 
function U (a, b, z) required for functional matching. One should use instead the form in 
Eq. (p.S|), which introduces logarithmic corrections to the power-law form of the solutions 
in regions I and II. These corrections then invalidate the matching condition in Eq. ( |3.'j| ) 
and in fact preclude the use of the functional matching procedure in the vicinity of the 
pole. Although one cannot a priori determine the values of /iM for which these logarithmic 
corrections are significant, the absence of peak- and pole-like features in both the numerical 
and point matching results, which do not depend on approximate expressions for U (a, 6, z), 
shows that the functional matching procedure does not give the correct values for the 
superradiant growth rate in a significant range of scalar masses. This is particularly evident 
for higher multipoles, as illustrated in Figure 4. 

Although the near and far regions are constructed in different ways in the present anal- 
ysis and in that of [25], these problems are inherent to the functional matching technique 
and must be taken into account in all applications of this procedure. Furthermore, the 
observed peak-like structures do not coincide exactly with those obtained in [25], occurring 
in the present analysis for smaller values of the scalar mass. Given that distinct anal- 
yses exhibit similar features at different points in parameter space, this further suggests 
that both peaks are simply the product of unsuitable approximations and should not be 
interpreted as physical enhancements of the superradiant growth rate. 

An extension of the point matching technique for a < M, although outside the scope 
of this work, would be useful for a better comparison with the numerical growth rates 
obtained in [13] close to the endpoint of the superradiant spectrum, as the continued- 
fraction method used in the latter work is designed for non-extremal black holes. Also, 
realistic black holes can only be spun up up to a ~ 0.998M, as estimated by Thorne in [32] 
taking into account the effects of radiation from the surrounding accretion disks. The 
extremal black hole bomb is nevertheless a very good approximation to the latter case, 
only differing significantly close to the endpoint of the spectrum, where the growth rate 
becomes very small. 

In spite of its limitations, the simple point matching method eliminates most of the 
shortcomings of other analytical computations of the superradiant spectrum and gives the 
best agreement with the numerical continued-fraction method. Despite its semi-analytical 
nature, as it requires solving the non-polynomial matching condition Eq. (|4.2| ), it is con- 
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siderably simpler than a full numerical solution and gives a good theoretical insight on 
the form of the bound state radial functions. This is crucial for several phenomenological 
purposes, in particular for determining the backreaction of the axionic clouds formed via 
superradiant emission around astrophysical black holes. This may, for example, affect the 
gravity wave signal produced by inspir ailing companions, which may lie within the reach 
of future observatories such as LISA^ [28]. 

The superradiant instability is extremely important in understanding the stability of 
rapidly-rotating black holes, at the same time providing a unique astrophysical probe of 
high-energy particle physics beyond the Standard Model, and we hope that this work 
motivates further exploitations of the rich phenomenological aspects of this mechanism. 
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